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Abstract 

Nuclear theory today aims for a comprehensive theoretical framework that can describe all nuclei. 
I discuss recent progress in this pursuit and the associated challenges as we move forward. 
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1. The journey toward a comprehensive theory of nuclei 

Nuclei comprise 99.9% of all baryonic matter in the Universe and are the fuel that 
burns in stars [T]. The rather complex nature of the nuclear forces among protons and 
neutrons generates a broad range and diversity in the nuclear phenomena that we ob- 
serve. Experiments indicate that developing a comprehensive description of all nuclei and 
their reactions requires theoretical and experimental investigations of rare isotopes with 
unusual neutron-to-proton ratios that are very different from their stable counterparts. 
These rare nuclei are difficult to produce and study experimentally since they can have 
extremely short lifetimes. The goal of a comprehensive description and reliable modeling 
of nuclei represents one of the great intellectual opportunities for nuclear physics in the 
coming years. 

Key scientific themes being addressed by the physics of nuclei and nuclear astrophysics 
communities are captured by five overarching questions: 

- What is the nature of the nuclear force that binds protons and neutrons into stable 
nuclei and rare isotopes? 

- What is the origin of simple patterns in complex nuclei? 

- What is the nature of neutron stars and dense nuclear matter? 
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- What is the origin of the elements in the cosmos? 

- What are the nuclear reactions that drive stars and stellar explosions? 

These questions align well with the drivers of rare isotope science as recently summarized 
in a report of the U. S. National Academies of Science Rare Isotope Scientific Assessment 
Committee (RISAC) [2]- One primary aspect of the first and second question concerns 
testing the predictive power of models by extending experiments to new regions of mass 
and proton-to-neutron ratio and identifying new phenomena that will challenge existing 
many-body theory. 

One of those challenges concerns the concept of nuclear shell structure and magic 
numbers. Understanding the origin of the classic, stable, nuclear magic numbers earned 
Goeppert-Mayer and Jensen the 1963 Nobel Prize. At the time, there were no data on 
very neutron-rich nuclei; however, theoretical predictions and experimental discoveries 
in the last decade indicate that nuclear shell structure is a rather fluid concept. For 
example, experimental data indicate that the magic numbers at N=20 and 28 fade away 
with neutron number, and the new magic numbers at N=14, 16, and 32 seem to appear 
(see, for example Ref. [3]). Nuclei far from stability also exhibit unusual properties as 
compared to their stable cousins. For example, the radial extension of the two-neutron 
halo in ^^Li (where ^Li acts as a core) is the same as that of ^°^Pb. 

We can conclude from just these two examples that nuclear structure is changing 
in the exotic environment. Since nuclei are self-bound, they generate their own mean 
field which determines how shell orbitals are filled. The experimental evidence indicates 
the way orbits are being filled is changing in neutron-rich nuclei. Furthermore, many- 
body correlations such as pairing superfluidity become crucial when the binding energy 
becomes small. These correlations also affect the open-shell character of nuclei, since the 
continuum of scattering states lies very close to the bound state. 

One of the interesting features of nuclear physics today is the overlaps and synergy 
between the physics of nuclei and nuclear astrophysics. The changes in nuclear structure 
occurring away from the valley of beta-stability do affect nucleosynthesis pathways. We 
have initial evidence of this on the measurements of the lifetime of ^*Ni [4] and its impact 
on r-process nucleosynthesis. Other examples involve the impact of electron capture on 
the core collapse mechanism of type-II supernovae :5! and the challenge of obtaining a 
nuclear equation of state that can quantitatively be used in predictions of neutron star 
properties 

A robust theoretical capability is required if we want to understand stable and exotic 
nuclei that are the core of matter and the fuel of stars. I will devote the remainder of this 
Proceedings to a discussion of some of the issues that remain in building this capability. 

2. Progress and challenges in the theory of nuclei 

The nuclear many-body problem spans nuclei from A=2 (the deuteron) to the super- 
heavy region, and different theoretical techniques are pursued in different regions of the 
chart of nuclei. This is depicted in Fig. 1 where the current reach of various theoreti- 
cal methods is interposed in the chart of nuclei. The challenges for theory can be viewed 
broadly as: determining the nuclear interaction, calculating ab initio nuclei, including the 
continuum, and finding an appropriate energy density functional for a nuclear Density 
Functional Theory approach. I will briefly discuss these areas in the following subsections. 
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Fig. 1. (Color online) The theoretical methods and computational techniques used to solve the nuclear 
many-body problem. The red vertical and horizontal lines show the magic numbers, reflecting regions 
where nuclei are expected to be more tightly bound and have longer half-lives. The anticipated path 
of the astrophysical r-process responsible for nucleosynthesis of heavy elements is also shown (purple 
line). The thick dotted lines indicate domains of major theoretical approaches to the nuclear many-body 
problem. For the lightest nuclei, ab initio calculations (Green's function Monte Carlo, no-core shell model, 
coupled cluster method), based on the bare nucleon-nucleon interaction, are possible (red). Medium- mass 
nuclei can be treated by configuration interaction techniques [interacting shell model (green)]. For heavy 
nuclei, the density functional theory based on self-consistent /mean- field theory (blue) is the tool of 
choice. By investigating the intersections between these theoretical strategics, one aims to develop a 
unified description of the nucleus. 

drawing from some of the work in which I have been involved. 

2.1. Determining nuclear interactions 

The forces among nucleons determine how nuclei behave. These forces arise from QCD. 
While there is some effort to extract directly from lattice calculations [7j the nuclear 
scattering, there have also been important developments in the application of a chiral 
Effective Field Theory (EFT) to the nuclear problem [8l9j . An EFT was developed for 
nuclei from a chiral Lagrangian that maintains all the symmetries of QCD and treats as 
fundamental particles the pions and nucleons. A power-counting scheme was developed 
so that the expansion could be carried out in a systematic way using as a small parameter 
of expansion the momentum transfer in nucleon-nucleon scattering divided by the QCD 
energy scale (1 GeV). This expansion has been carried out to fourth-order (known as 
N'^LO) [TU]. The resulting nucleon-nucleon potential fits all scattering phase shifts up to 
the cut-off used in the EFT. 

An interesting feature of the EFT is that three-nucleon interactions arise naturally at 
the third order (N^LO). Chiral three-body forces have been known for a long time, but 
their systematic computation was not known until recently. Nuclear spectroscopy using 
three-body interactions from the EFT expansion improve when compared to calculations 
with two-body forces alone [TT| . This is also to be expected since it has long been known 
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that two-body forces alone do not reproduce the correct binding in ^He. 
Thus, our starting point is a Hamiltonian of the form 

H = t- i„„ + 1^2 + 1^3 , (1) 

where t is the kinetic energy operator and tcm is the center-of-mass operator (so that 
the kinetic energy maintains translational invariance) and the V2 and V3 represent the 
two-body and three-body nuclear potentials. The three-body interaction can actually be 
written in normal-ordered form as 

1/3 = i ^{ijk\\ijk) + ]^'^{ijp\\ijq){alaq} + {ipq\\irs){ala\asar} + h , (2) 

ijk ^jPQ ipqrs 

where {aP"f\\Ser]) represents a matrix element of the three-body interaction, i, j, k denote 
the hole states, a, b, c denote the particle states, p, g, r denote all states, and denotes 
the residual three-body Hamiltonian 

^3 = ^ X! (wl|s^w){apajaj;a„ataj . (3) 

pqrstu 

This notation will be useful in what follows. 

V2 is often modeled with a very repulsive short distance interaction. In the case of 
Avl8 T2], the central repulsion at r = is approximately 3 GeV; in N^LO, it is 500 
MeV [10 . In order to use these interactions efficiently in a framework where one expands 
the many-body wave function in a basis, it is desirable to renormalize these interactions. 
This can be done in several ways, including the G-matrix approach [13] and similarity 
transform methods ,14J. 

More recently, work in the direction of renormalization group (RG) [15|16j techniques 
has been used to evolve the nucleon-nucleon interactions from high momentum to low 
momentum while preserving the nucleon-nucleon phase shifts up a given momentum. 
These interactions can then be projected onto the appropriate basis states for computa- 
tion. These approaches offer promise with a price: the resulting two-body interactions do 
not recover the binding energy that would be obtained from a pure two-body calculation, 
and one must indeed include a three-body force in the calculations. Furthermore, the RG 
interactions tend toward overbinding both in nuclei |17) and in nuclear matter [18j . The 
cure to this problem is the introduction of three-body forces which naturally occur in the 
RG approaches. Another possibility is to adjust the three-body force that comes from 
EFT to obtain the experimental binding of "^He [19] . 

2.2. Calculating ah initio Nuclei 

For light nuclei with mass numbers A < 12, both Green's function Monte Carlo meth- 
ods [20 and no-core shell-model calculations |21lllj using a basis-state expansion pro- 
vide almost converged benchmarks for selected two- and three-body Hamiltonians. The 
agreement with experimental data for many light nuclei is quite reasonable in these cal- 
culations. GFMC techniques do not need to utilize renormalized interactions, but they 
can only use local potentials, and the main group involved uses one two-body potential 
(Avl8, [12j ) and an appropriate three-body potential fitted to experimental data in light 
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nuclei 22]. The NCSM approach can use any potential, but in most cases uses renormal- 
ized ones in order to obtain converging calculations. Both methods exhibit exponential 
scaling with the number of particles (to reach the same accuracy as in "^He with the same 
methods) . Recently, we have seen a movement of these methods into the realm of light 
nuclear reactions. For example, the GFMC approach was used to describe n + a scat- 
tering |23j while the NCSM was used to calculate the ^Be(p, 7)^B S-factor [24]. Recent 
progress has also been made in ab initio inelastic four-body scattering problems pS] , 

Coupled-cluster theory [26l27l28l29j represents another approach to ab initio calcu- 
lations of nuclei. Along with collaborators, including computational chemists, we have 
developed the theory and methods |30|31|32|33|17j to the point where we can calcu- 
late nuclear properties in medium-mass nuclei. Nuclear coupled-cluster theory exhibits 
polynomial scaling with the size of the basis state and number of particles, and is size- 
extensive. This very important property means that no unlinked diagrams enter into the 
theory. 

Coupled-cluster theory starts with the simple assertion that the correlated ground- 
state wave function can be described by applying an exponentiated correlation operator 
to an uncorrelated Slater determinant that naively describes the nucleus: 

I = exp (T) I $ ) . (4) 

The correlation operator T induces various np-nh (particlc-hole) correlations (up to the 
number of nuclcons. A, in the nucleus): 

T = Ti+T2 + ---Ta. (5) 

The energy of the ground state is given by 

E={^\ exp(-T)ff exp(r) | $) = ($ | | $) = ($ | {Hexp(T))^ \ $) , (6) 

where the subscript c means that only connected diagrams enter into the expressions. 
The T„ operators take the form (as is the case here for the two-body excitation operator) 

^2 = E tfj^aaUj^^ (7) 

ab^ij 

where tfj are the 2p-2h correlation amplitudes. 

One must compute the correlation amplitudes through equations that left-project H 
onto the space of excited Slater determinants. For example, if the theory only contains 
Ti and T2 operators, then the equations to solve for the coefficients tf and tf^ are 

= I I $) 

= ($^/|iIl$), (8) 
and so on. 

Eqn. ^ is exact; however, the power of the coupled-cluster theory is its highly accurate 
computation of the energy even when one limits to lower order the number of T„ operators 
in the theory. For example, at the Ti and T2 levels (coupled clusters in singles and doubles, 
or CCSD), approximately 90% of the correlation is obtained, while with approximate 
triples corrections (denoted CCSD(T)), nearly all of the correlation energy is obtained 
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-'^CCSD{T) 


-29.2 


-148.2 


-502.9 


exact (FY) 


-29.19(5) 







Table 1 

Reference vacuum energies, Eg, CCSD and CCSD(T) correlation energies, Ai?ccSD a-nd Ai?cQSD{T)i 
and binding energies i?ccSD{T) ''He, and '^'^Ca. The vacuum energies, Eq, are for hu> = 14 MeV 
in the case of "He and = 22 MeV for ^''O and *°Ca. The CCSD and CCSD(T) energies are the 
extrapolated infinite model space results. The exact Faddeev-Yakubovsky result is from Ref. |19) 

|31ll7j . This contrasts sharply to truncated shell-model calculations where the truncation 
introduces unlinked diagrams causing growing errors with particle number [34j . 

A certain amount of benchmarking of methods is necessary in order to understand 
their regions of validity. For coupled-cluster theory, this has recently been performed in 
the triton and '^He. We also performed extremely large basis-set calculations to which 
other methods can compare when they reach the technical capability to do so [17]. We 
summarize our coupled-cluster results for the binding energies of ''He, ^^O, and ''"Ca 
in Table [1] which gives the extrapolated correlation energies Ai?ccSD and Ai<^ccSD(T) j 
relative to the uncorrelated energy Eq = \ H \ $}. We find that for ^He, ^^O, and 
^°Ca, the triples corrections are a factor of w 0.015, 0.066, and 0.081 smaller than the 
CCSD correlation energies. From this, we estimate the missing correlation energy from 
quadruples, pentuplets, and so on to be of the order of 1 MeV for "'"Ca. We note that 
^^O is overbound by about 20 MeV, and "'"Ca by about 150 MeV when compared to the 
experimental binding energies. This is not surprising and points to the importance of 3NF 
for nuclear structure calculations and to its apparently large role in the RG prescriptions. 

The three-body force is a key component of modern nuclear Hamiltonians. Other meth- 
ods have been adapted to include three-body interactions already several years ago. We 
recently developed coupled-cluster theory for three-body forces [S^- This represents a 
major improvement to the scientific reach of modern coupled-cluster calculations for nu- 
clei. An exponential extrapolation of the (approximate) CCSD(T) minima to an infinite 
model space is shown in Fig. 2 and yields E^o — —28.24 MeV. This is in excellent agree- 
ment with the exact Faddeev-Yakubovsky result E = —28.20(5) MeV. In our largest 
model space at the minimum hiu = 17 MeV, the ground-state expectation values for 
the center-of-mass Hamiltonian is (i?cm) ~ 20 keV while the expectation value for the 
angular momentum is zero for a closed-shell nucleus by construction. 

It is interesting to analyze the different contributions AE to the binding energy E. 
The individual contributions are given in Fig. 2 for a model space of iV = 4 oscilla- 
tor shells and Hco = 20 MeV. The main contribution stems from the low-momentum 
NN interaction. The contributions from 3NFs account only for about 10% of the total 
binding energy. This is consistent with the chiral EFT power-counting estimate (Vsn) ~ 
{Q/K^Y (^lowkfe) ~ 0.1 (Mowfe) [Tpl. The second-, third-, and fourth-largest contribution 
are due to the first, second, and third term on the right-hand side of Eq. ([2]). These are 
the density-dependent zero-, one-, and two-body terms, which resulted from the normal 
ordering of the three-body Hamiltonian in coupled-cluster theory. The contributions from 
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Fig. 2. Left: (Color online) Data points: CCSD(T) results (taken at the hui minima) for the binding 
energy of ^He with 3NFs as a function of the number of oscillator shells. Dashed lines: Exponential fit to 
the data and asymptote of the fit. Full line: Exact result. Right: Relative contributions \AE/E\ to the 
binding energy of ^He at the CCSD level. The different points denote the contributions from (1) low-mo- 
mentum NN interactions, (2) the vacuum expectation value of the 3NF, (3) the normal-ordered one-body 
Hamiltonian due to the 3NF, (4) the normal-ordered two-body Hamiltonian due to the 3NF, and (5) 
the residual 3NFs. The dotted line estimates the corrections due to omitted three-particle— three-hole 
clusters. 

the residual three-body HamUtonian, Eq. are very small and are represented by the 
last point in the right panel of Fig. 2. The residual 3NF contributes to the energy directly 
and indirectly through a modification of the Ti and T2 cluster amplitudes. Apparently, 
both contributions are very small. This is a promising result, and time will tell whether 
it holds for heavier nuclei. 



2.3. Presence of the continuum 

Exotic phenomena emerge in weakly bound and resonant many-body quantum sys- 
tems. These phenomena include ground states that are embedded in the continuum, 
melting and reorganizing of shell structures, extreme matter clusters and halo densities. 
A theoretical description of weakly bound and unbound quantum many-body systems 
represents a challenging undertaking. The weakly bound character of these systems means 
that they should be treated as open quantum systems where coupling with the scattering 
continuum can take place. Furthermore, theoretical treatments that consider continuum 
basis states can describe resonant widths (lifetimes). 

Recent work with Gamow states employed in Hamiltonian diagonalization methods 
[35136137] have shown that these basis states correctly depict properties associated with 
open quantum systems. This Berggren basis is composed of bound, resonant, and (contin- 
uum) scattering single-particle states |38| . This basis significantly improves and facilitates 
the description of loosely bound systems and is essential in the description of unbound 
systems. In addition, several groups have worked on alternative methods, such as the 
continuum shell model |39l40j and the recently developed shell model embedded in the 
continuum 41 . 

We recently developed a complex version of the coupled-cluster method that can utilize 
the Gamow-basis to calculate widths of states in the He isotopic chain in an ab initio 
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Fig. 3. (Color online) CCSD calculation of the ''"^"He ground states with the low-momentum N'^LO 
nucleon-nucleon interaction for increasing number of partial waves. The energies E are given in MeV for 
both real and imaginary parts. Experimental data are from Ref. 42'. Our calculated width of ^"He is « 
0.002 MeV. 



framework [32]. We utilized a Viowfc interaction from N'^LO at a cut-off of 1.9 fm^^ in a 
Berggren basis to perform the calculations. These novel calculations, the results of which 
are shown in Fig. 3, while by no means perfect, utilizing only a two-body interaction, do 
indicate the power of moving beyond the shell model. 



2.4. Density functional theory for nuclei 

Nuclear Density Functional Theory (nDFT) is being used and improved to describe 
heavier nuclei. The basic idea is to expand the energy density functional to second order 
(gradient terms) in the density and to calculate with this functional all nuclear proper- 
ties. One must determine the coupling constants that characterize the functional from 
experimental data. Odd-A and odd-odd nuclei will be key to any new parameterizations 
as they allow one to interrogate time-odd terms in the energy density functionals. 

The nuclear DFT efforts have been quite successful in describing a wide variety of 
nuclear data with very good precision across the nuclear chart. The various parameteri- 
zations usually work quite well in regions where nuclear masses and other properties are 
experimentally determined, but extrapolations into very neutron-rich nuclei have been 
problematic. The next-generation experimental facilities should enable theorists to ob- 
tain a functional parameterization that will describe bulk properties of all nuclei. Various 
nuclear data along long isotopic and isotonic chains are needed to constrain the isovector 
part of the energy functional. More specifically, one needs (difference of) masses and 
measures of collectivity and of the shell evolution in unknown regions, where predictions 
of currently used functionals disagree. Initial attempts are represented in Fig. 4 where 
the two-neutron separation energies were calculated from an even-even mass table for the 
Sly4 Skyrme interaction with pairing included and partial number projection |43;. Other 
groups are also investigating nuclear mass systematics in large mass-table calculations 
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Neutron Njmb&r 



Fig. 4. (Color online) An example of a large-scale systematic density functional theory (DFT) calculation 
for complex nuclei: results of the deformed DFT calculations of two-neutron separation energies for 1553 
particle-bound even-even nuclei with Z < 108 and N < 188 (from M. Stoitsov). 



and have recently extended these calculations to include a study of fission barriers [33] . 
Other data that could be included in determining the coupling constants of the energy 
density functional, such as large deformations (at low and high angular momentum) and 
multipole strength distributions in neutron-rich nuclei, will also be extremely valuable. 
The quest for a universal energy density function represents one of the key challenges of 
a U.S. Department of Energy SciDAC project. 



3. The next 10 years 

I believe that the next 10 years hold great promise for the physics of nuclei and for 
moving toward the ultimate solution of the nuclear problem: to obtain a comprehensive 
understanding of all nuclei and their interactions. On the experimental side, new facilities 
at RIKEN, GSI, and GANIL, and the future FRIB in the U.S. will make the discoveries 
that will challenge our theories and our understanding of many-body phenomena. 

From a theoretical point of view, I believe the key to reaching the goal is the pursuit of 
many-body theory and its connections to nuclear DFT. We will see in the next few years 
the development of nuclear theory such that robust and predictive ab initio calculations 
can be accomplished in medium-mass nuclei. The effort to establish an improved density 
functional should allow for breakthrough science in heavy nuclei and in nuclei along the 
various nucleosynthesis paths. The strides being made to connect ab initio calculations 
in light- to medium-mass systems to nuclear DFT will enable us to finally rest the energy 
density functional on a firmer footing. The links from QCD to the underlying interaction 
will allow us to close the loop and reliably compute nuclei from the ground up. 

All of the efforts I described in the preceding section ~ GFMC, NCSM, coupled-cluster 
efforts, and the quest to find the universal nuclear energy density functional - require 
petascale computational capability. We have witnessed in the last five years tremendous 
growth in computing capability and many in our community continue to embrace these 
marvelous tools. The ingenuity of researchers to utilize the largest of these machines 
undoubtedly will propel the field forward in our quest to understand nuclei. 
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